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> ■ ABSTRACT 

Aims. We study degeneracies between cosmological parameters and measurement errors from cosmic shear surveys. We simulate 
realistic survey topologies with non-uniform sky coverage, and quantify the effect of survey geometry, depth and noise from 
I intrinsic galaxy ellipticities on the parameter errors. This analysis allows us to optimise the survey geometry. 
Methods. We carry out a principal component analysis of the Fisher information matrix to assess the accuracy with which 
linear combinations of parameters can be determined. Using the shear two-point correlation functions and the aperture mass 
dispersion, which can directly be measured from the shear maps, we study various degeneracy directions in a multi-dimensional 

■ parameter space spanned by f2m, ^a, erg, the shape parameter F, the spectral index ris, along with parameters that specify the 
Oh' distribution of source galaxies. 

^ , Results. A principal component analysis is an effective tool to probe the extent and dimensionality of the error ellipsoid. If 

5^ ' only three parameters are to be obtained from weak lensing data, a single principal component is dominant and contains 
' all information about the main parameter degeneracies and their errors. For four or more free parameters, the first two 
I principal components dominate the parameter errors. The degeneracy directions are insensitive against variations in the noise 

1^ or survey geometry. The variance of the dominant principal component of the Fisher matrix, however, scales with the noise. 
. ^ [ Further, it shows a minimum for survey strategies which have small cosmic variance and measure the shear correlation up to 
• several degrees. This minimum is less pronounced if external priors are added, rendering the optimisation less effective. The 

$_( minimisation of the Fisher error ellipsoid can lead to slightly different results than the principal component analysis. 
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1. Introduction 

Recent observations by the Wilkinson Microwave 
Anisotropy Probe (WMAP) mission confirmed the 
standard cosmological model with a very high degree of 
accuracy (Spergel et al. 2003). In particular, these obser- 
vations confirmed that the universe is spatially flat and 
dominated by dark energy and dark matter. Regarding 
the initial power spectrum of scalar perturbations, the 
predictions of the simplest inflationary models were 
strengthened, i.e. the near scale-invariance, adiabaticity 
and Gaussianity of the initial density perturbations. 
However, certain outstanding issues remain to be solved. 
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such as the running of spectral index as which can be 
addressed in more detail with additional data from galaxy 
surveys such as SDSS (York et al. 2000), 2dF (CoUess et 
al. 2001) and the Lyman-a forest (see Seljak et al. 2003 
and references therein). 

Weak lensing surveys are expected to make impor- 
tant and complementary contributions to high-precision 
measurements of cosmological parameters. Contaldi et al. 
(2003) used the Red Cluster Sequence (RCS) to show 
that the r^m-CTs degeneracy direction is nearly orthogo- 
nal to the one from CMB measurements, making weak 
lensing particularly suitable for combined analyses (van 
Waerbeke et al. 2002). Ishak et al. (2003) argued that a 
joint CMB-cosmic shear survey provides an optimal data 
set for constraining the amplitude and running of spec- 
tral index which helps to probe various inflationary mod- 
els. Tereno et al. (2004) studied cosmological forecasts 
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for joint CMB and weak lensing data. Clearly, the po- 
tential of weak lensing surveys (Mellier 1999; Bartelmann 
& Schneider 2001; Refregier 2003; van Waerbeke & Mellier 
2003; Schneider 2005) as a cosmological probe is now well 
estabhshed (Contaldi et al. 2003; Hu & Tegmark 1999). 
In the last few years there have been many studies which 
have detected cosmic shear in random patches of the sky 
(Brown et al. 2003; Bacon et al. 2003; Bacon, Refregier 
& Ellis 2000; Hamana et al. 2003; Hammerle et al. 2002; 
Hoekstra et al. 2002a; Hoekstra, Yee & Gladders 2002a; 
Jarvis et al. 2002; Kaiser, Wilson & Luppino 2000; Maoh 
et al. 2001; Refregier, Rhodes, & Groth 2002; Rhodes, 
Refregier & Groth 2001; van Waerbeke et al. 2000, 2001, 
2002; Wittman et al. 2000). While early studies were pri- 
marily concerned with the detection of the weak lensing 
signal, present generations of weak lensing observations 
are putting constraints on cosmological parameters, in 
particular the matter density parameter f2„i and the power 
spectrum normalisation erg. 

Inspired by the success of these surveys, there are many 
other ongoing, planned and proposed weak lensing sur- 
veys which are currently in progress, including the Deep 
Lens Survey (Wittman et al. 2002), the Canada- France- 
Hawaii Telescope Legacy Survey (Hoekstra et al. 2005; 
Semboloni et al. 2005), the Panoramic Survey Telescope 
and Rapid Response System, the Supernova Acceleration 
Probe (Massey et al. 2003), the NOAO Deep Wide-Field 
Survey (Groch et al. 2002) and the Large Synoptic Survey 
Telescope (Tyson et al. 2002). Future cosmic shear sur- 
veys will be able to probe much larger scales in the linear 
regime which will provide more stringent bounds on cos- 
mological parameters such as the equation of sate of dark 
energy and its time variations. 

In a recent work (Kilbinger & Schneider 2004, hereafter 
KS04), the impact of the survey design on cosmological 
parameter constraints was analysed using a likelihood and 
Fisher matrix analysis, extending previous studies based 
on the assumption of uniform sky coverage (Schneider et 
al. 2002a). Earlier work in this direction by Kaiser (1998) 
considered a singe 3° x 3° -field and studied the effect of 
sparse sampling and intrinsic ellipticity dispersion. The 
motivation for the present work remains the same, al- 
though we concentrate on the eigenvalues of the Fisher 
matrix. Therefore, in contrast to KS04 where the la errors 
on individual parameters have been used to optimise the 
survey geometry, we consider all parameter combinations 
corresponding to the eigenvectors of the Fisher matrix. 
We study how noise due to the intrinsic ellipticity disper- 
sion of galaxies a^, the number density of galaxies rigai, 
the survey depth and marginalisation affects the determi- 
nation of the parameter combinations for various survey 
strategies. 

Cosmic shear is sensitive to a large number of cosmo- 
logical parameters. However, the dependency on these pa- 
rameters is partially degenerate (although these degenera- 
cies can be broken by the use of external data sets such as 
CMB, galaxy surveys and Lyman-a surveys). A principal 
component analysis (PCA) can be used as an efficient tool 



to identify the degeneracy directions and linear combi- 
nations of cosmological parameters, rank-ordered accord- 
ing to the accuracy with which they can be determined 
from a given survey set-up. Indeed, in recent years there 
has been a renewed interest in applying principal com- 
ponent analysis techniques to various cosmological data 
sets, a technique pioneered by Efstathiou & Bond (1999). 
This method can reveal the detailed statistical structure 
of cosmological parameter space which is lacking in an 
one-dimensional confidence level presentation. Efstathiou 
(2002) studied PCA in the context of the tensor degener- 
acy in CMB. For a recent work see Rocha et al. (2004), 
where the possibility of measurement of the fine-structure 
constant a has been explored in the context of CMB data 
with analysis based on Fisher matrix and PCA. Hu &; 
Keeton (2002) applied this technique to map the density 
distribution along the radial direction from weak lensing 
surveys. Jarvis & Jain (2004) used PCA to correct for the 
point spread function (PSF) variation in weak lensing sur- 
veys. In the context of SN la observations to constrain the 
dark energy equation of state, Huterer & Starkman (2003) 
and Huterer & Cooray (2004) employed PCA and its vari- 
ants (see Wang & Tegmark 2005; Crittenden & Pogosian 
2005 for more recent results). Tegmark et al. (1998) used 
PCA for decorrelating the power spectrum of galaxies. 
This idea was initially proposed by Hamilton (1997) and 
further discussed in the context of galaxy surveys by 
Hamilton & Tegmark (2000). 

This paper is organised as follows: in Sect. 13 a very 
brief overview of our notations is provided; in particular, 
we introduce how the covariance matrix and the Fisher 
matrix is constructed for a given estimator and a given sur- 
vey strategy. This section also outlines the basics of prin- 
cipal components analysis. In the next section (Sect.O we 
provide the details of survey geometries and the numerical 
results of the PCA considering three survey set-ups. For a 
small number of cosmological parameters, (< 4) we con- 
sider various survey strategies and try to optimise those. 
For a larger set of parameters, we consider a ten times 
larger shear survey. Effects of various noise sources on the 
principal components are investigated. Sectionals left for 
discussions and future prospects. 



2. Notation and formalism 

2.1. Second-order shear statistics 

In our numerical studies presented here, we use the two- 
point correlation functions of shear ^j- and the aperture 
mass dispersion (Mfp) to predict constraints on cosmolog- 
ical parameters. Both these statistics depend linearly on 
the convergence power spectrum (Kaiser 1992; Kaiser 
1995; Schneider 1996; Schneider et al. 1998) 
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where Ji, is the first-kind Bessel function of order v. 

Estimators of these statistics and their covariances are 
defined in Schneider et al. (2002a). We use the Monte- 
Carlo-hke method from KS04 to integrate the analytical 
expressions of the covariances, which are exact in case of 
a Gaussian shear field. Our result is expected to underes- 
timate the covariance due to non-Gaussian contributions 
on scales between ~ 1 and 10 arc minutes. 

2.2. Fiducial cosmological model 

We calculate the convergence power spectrum and the 
shear estimators using the non-linear fitting formulae of 
Peacock & Dodds (1996). Our cosmological model has 
seven free parameters: These are the five cosmological pa- 
rameters firm ^A, the power spectrum normalisation ag, 
the spectral index of the initial scalar fluctuations rig and 
r, which determines the shape of the power spectrum. The 
fiducial model is assumed to be a flat ACDM cosmology 
with ilra = 0.3, erg = 1, r = 0.21 and Us = 1. The two 
parameters zo and (3 characterise the redshift distribution 
of background galaxies (Brainerd et al. 1996), 

with fiducial values zq — I and (3 = 1.5. 

2.3. Principal components analysis of the Fisher matrix 

We use the expression for the Fisher matrix (see Tegmark, 
Taylor & Heavens 1997 for a review) in the case of 
Gaussian errors and parameter-independent covariance. 




where Xk is either C+(^'fc)> i-i^k), (^ap(^'fc)) or an entry 
of the combined correlation function ^tot = i^+T^-)- C 
denotes the covariance matrix of the estimator of the cor- 
responding shear statistics, = (Oi, . . . 0„) is the vector 
of cosmological parameters. The inverse of the Fisher ma- 
trix is the covariance of the parameter vector at the point 
of maximum likelihood, 

= (A0A0*) = (00*) - (0)(0*). (4) 

The standard deviation of the z**^ parameter obtained from 
the Fisher matrix, AG, = ((ef)-(e,)2)i/2 = [(F-i),,]i/2^ 
is called the minimum variance bound (MVB). According 
to the Cramer-Rao inequality, the variance of any unbi- 
ased estimator is always larger or equal to the MVB. 

Any real matrix W is called a decorrelation matrix if 
it satisfies 

F = W*AW, (5) 



where A is a diagonal matrix (Hamilton & Tegmark 2000). 
The quantities $ = W0 are uncorrelated because their 
covariance matrix is diagonal, 

(A*A*') = W(A0A0')W' = A-^ (6) 

By multiplying W with the square root of the diagonal 
matrix A, the quantities $ can be scaled to unit variance 
without loss of generality. In this case, © is written as 

F = W'W, (7) 

where W = A^/^ W. Note that the choice of W is not 
unique. If some matrix W satisfies {TJ, the same is true 
for any orthogonal rotation OW with O G SO(n) and 
therefore, there are infinitely many decorrelation matrices 
satisfying (jZj). 

If W is an orthogonal matrix, its rows are the eigen- 
vectors Pi of F and A = diag(Ai) is the diagonal matrix of 
the corresponding eigenvalues. In that case, © is called a 
principal component decomposition. We assume the eigen- 
values to be in descending order. 

The eigenvectors or principal components of F deter- 
mine the principal axes of the n-dimensional error ellipsoid 
in parameter space. The eigenvectors represent orthogonal 
linear combinations of the physical (cosmological) parame- 
ters that can be determined independently from the data. 
The more these vectors are aligned with the parameter 
axes, the less are the degeneracies between those parame- 
ters. The accuracy with which these linear parameter com- 
binations can be determined is quantified by the variance 
Oi EE cr(p,) = A$, = A^^/^ = A^^/^ Thus, a princi- 
pal component decomposition of the Fisher matrix gives 
us information about which (linear) parameter combina- 
tions can be determined with what accuracy from a given 
data set. Since the eigenvalues are in descending order, 
the first eigenvector pi having the smallest variance cor- 
responds to the best constrained parameter combination. 
The last eigenvector p„ is the direction with the largest 
uncertainty. 

From 0, one can calculate the MVB from the eigen- 
vectors and eigenvalues of F, 

^Q,^\Y.wlK'\ ■ (8) 

3. Numerical results 

3.1. Survey strategies 

We simulate shear surveys consisting of P circular, uncor- 
related patches of radius R on the sky, each in which N 
individual fields of view of size 13' x 13' arc distributed 
randomly but non-overlapping. The total number of fields 
of view is 71 = PN — 300, corresponding to a total survey 
area of ^ = 14.1 square degree. Different surveys with 
N — 10, 20, 30, 50, and 60 are considered, correspond- 
ing to geometries with P = 30, 15, 10, 6 and 5 patches, 
respectively (see KS04). 
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We denote these survey strategies with {N,R), e.g. 
(50, 100') corresponds to a survey with TV = 50, i? = 100' 
and P = 6. Further, we consider a survey consisting of 300 
uncorrelated hues of sight a 13' x 13', which are randomly 
distributed on the sky. This survey, denoted by 300 -13'^, 
has smaller cosmic variance than any of the patch strate- 
gies, but does not sample intermediate and large angular 
scales. 

If not indicated otherwise, the number density of 
source galaxies is ngai = 30arcmin~^. This number 
density of high-redshift galaxies which are usable for 
weak lensing shape measurements can be achieved with 
high-quality ground-based imaging data on a 4 m-class 
telescope. The source galaxy ellipticity dispersion is 
(Te = 0.3, if not stated otherwise. For comparison, these 
quantities are varied to rigai — 20 and 40arcmin^^, and 
tie = 0.2, 0.4, respectively, to study the effect of noise 
sources on the principal components. 
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3.2. Eigenvalues of the Fisher matrix 

We consider the Fisher matrix F corresponding to all 
seven cosmological and redshift parameters. In this sec- 
tion, we demonstrate the influence of survey characteris- 
tics (other than the geometry) on the eigenvalues of F. 
The variance of the linear combination of parameters given 
by the i'^ eigenvector pi of F is cr(pi) = ^i^^^, as defined 
in Sect. 12.31 In Fig. ^ we show the effect of the number 
of background galaxies, rigai, and the intrinsic ellipticity 
dispersion, a^. In both cases, the noise variation causes a 
scaling of the variance. For the intrinsic ellipticity disper- 
sion, the scaling factor increases with i. The i = 1 vari- 
ance scales linearly with , whereas the mean dependence 
(averaged over all 7 eigenvalues) is quadratic in ■ In the 
case of ngai, however, the variance <Ti of all eigenvectors is 
scaled by a constant factor which is inversely proportional 
to rigai. The variance of the eigenvectors are steeper func- 
tions of the noise characteristics than the MVB (Kilbinger 
& Munshi 2005). Note however, that in this previous study 
the MVB was calculated for each parameter individually, 
without taking parameter correlations into account. 

If boundary effects due to the finite field of the survey 
are neglected, the covariance is anti-proportional to the 
observed survey area. Consequently, the variance (Ji scales 
as /si;y^^, where /sky is the fractional sky coverage of the 
survey. As an example, we compare the 300 • 13'^ survey 
with a survey consisting of 50 patches with = 60 and 
R = 140' (corresponding to ten (60, 140') surveys). We 
found a good agreement on the expected scaling of (7^ as a 
function of the survey area, although for the combinations 
which are worst constrained by the data, the dependence 
on /sky seems to be steeper (see Fig. EJ. 



Fig. 1. The variance a{pi) — associated with the 

principal components of the 7x7 Fisher matrix. The 
straight lines are fits to the data points. Lower left panel: 
<y{pi) as a function of the intrinsic ellipticity dispersion of 
galaxies (Jg = 0.2 and 0.4. Lower right: the variation of 
<7{pi) with the number density of galaxies rigai = 20 and 
40. Upper left: the variation of cr(p„) for three different 
estimators, ^_ and ^tot- Upper right: the variation of 
(j{Pi) for three different priors, see Sect. 13.111 The survey 
strategy is (30, 100') in all cases. 




Fig. 2. The variation of cr(pi) as a function of the eigen- 
vector number i, for ^tot (left panel) and (Map) (right 
panel). The redshift parameter f3 is fixed, and a Gaussian 
prior with variance s{ag,) = 0.1 is added to regularise the 
Fisher matrix. Two surveys with size 14.1 (solid lines) and 
141 square degree (dashed), respectively, are displayed. 



3.3. Eigenvectors of the Fisher matrix 

In Tables the eigenvectors of the Fisher matrix cor- 
responding to various combinations of parameters are 
shown, for the 2PCF and the aperture mass dispersion. 
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Table 5. Eigenvectors of the Fisher matrix corresponding to (fim, T, erg, rig) and the prior Oa = 0.7, for ^^^t 
in brackets) using the survey strategy (50, 100'). 





pi 


P2 


P3 


P4 






0.657 


-0.588 


0.282 


-0.376 


0.154 




(0.7Uo) 


(—0.533) 


(U.320) 


(—0.332)+ 


(0.158)+ 




(0.591)^ 


(—0.647)- 


Co 3291- 


(—0.349)- 


(n 290)~ 

\ yj < ^ \j J 


r 


0.474 


0.719 


-0.217 


-0.459 


0.182 




(0.406)+ 


(0.768)+ 


(-0.116)+ 


(-0.480)+ 


(0.198)+ 




(0.547)" 


(0.654)" 


(-0.205)" 


(-0.480)- 


(0.389)- 




0.542 


-0.062 


-0.474 


0.690 


0.279 




(0.548)+ 


(0.005)+ 


(-0.572)+ 


(0.611)+ 


(0.289)+ 




(0.533)- 


(-0.129)- 


(-0.522)- 


(0.653)- 


(0.538)- 


Us 


0.217 


0.363 


0.805 


0.415 


0.209 




(0.179)+ 


(0.352)+ 


(0.746)+ 


(0.535)+ 


(0.290)+ 




(0.257)" 


(0.367)" 


(0.760)" 


(0.470)- 


(0.414)- 




0.004 


0.013 


0.165 


0.389 






(0.005)+ 


(0.016)+ 


(0.255)+ 


(0.408)+ 






(0.006)- 


(0.020)- 


(0.221)- 


(0.804)- 





Table 6. Eigenvectors of the Fisher matrix corresponding to (Om, F, erg, Oa), for ^tot (^+)$- in brackets) using the 
survey strategy (50, 100'). 





pi 


P2 


P3 


Pi 


AQj 




0.676 


0.556 


0.082 


-0.475 


0.147 




(0.722)+ 


(0.501)+ 


(0.030)+ 


(-0.475)+ 


(0.172)+ 




(0.613)- 


(0.619)- 


(0.348)- 


(-0.343)- 


(0.213)- 


r 


0.479 


-0.773 


0.385 


-0.156 


0.094 




(0.408)+ 


(-0.808)+ 


(0.371)+ 


(-0.207)+ 


(0.112)+ 




(0.558)- 


(-0.718)- 


(0.400)- 


(0.110)- 


(0.156)- 


(T8 


0.555 


0.022 


-0.339 


0.759 


0.244 




(0.556)+ 


(-0.039)+ 


(-0.260)+ 


(0.788)+ 


(0.290)+ 




(0.549)" 


(0.088)" 


(-0.721)" 


(0.441)" 


(0.333)- 




-0.060 


0.303 


0.854 


0.417 


0.220 




(-0.032)+ 


(0.307)+ 


(0.891)+ 


(0.332)+ 


(0.232)+ 




(-0.092)" 


(0.304)" 


(0.445)" 


(0.837)" 


(0.448)- 


^-1/2 


0.005 


0.014 


0.210 


0.308 






(0.005)+ 


(0.017)+ 


(0.223)+ 


(0.362)+ 




(Ar^/^)- 


(0.006)- 


(0.021)- 


(0.364)- 


(0.499)- 





Table 1. Eigenvectors of the Fisher matrix corresponding 
to (ilm, r, as) and a flat Universe as prior, for (Mfp) using 

the survey strategy (50, 100') . ' = CTj is the variance of 
the i"^ eigenvector, Qj the MVB for the j**^ cosmological 
parameter. 





pi 


P2 


P3 


Ae^ 




0.649 


-0.533 


-0.542 


0.165 


r 


0.563 


0.816 


-0.128 


0.042 




0.511 


-0.222 


0.830 


0.252 




0.004 


0.020 


0.300 





Table 2. Eigenvectors of the Fisher matrix corresponding 
to (Jim, r, (Tg, Ug) and the prior f^A = 0.7, for (M^p) using 
the survey strategy (50, 100'). 





pi 


P2 


P3 


P4 


AG, 




0.563 


-0.671 


0.322 


-0.358 


0.254 


r 


0.575 


0.622 


-0.235 


-0.474 


0.333 


0-8 


0.522 


-0.161 


0.501 


0.670 


0.473 


Us 


0.281 


0.368 


0.767 


0.443 


0.329 




0.004 


0.015 


0.147 


0.698 





The best determined eigenvector pi is always orthogo- 
nal to the fim-cg degeneracy direction. The variance (t„ of 



the worst constrained principal component dominates the 
uncertainty of all eigenvectors. In the case of (fim, E), 
af = I/A3 constitutes more than 99% of the total un- 
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Table 3. Eigenvectors of the Fisher matrix corresponding 
to (fim, r, o-gj ^a) for {M'ip) using the survey strategy 
(50,100'). 





Pi 


P2 


P3 


P4 


Ae, 




0.587 


0.643 


0.352 


-0.341 


0.169 


r 


0.591 


-0.694 


0.392 


0.114 


0.112 


0-8 


0.542 


0.117 


-0.727 


0.404 


0.252 




-0.102 


0.299 


0.439 


0.840 


0.368 




0.004 


0.015 


0.258 


0.417 





Table 4. Eigenvectors of the Fisher matrix corresponding 
to (flm, r, as) and a flat Universe as prior, for ^tot (C+i 
in brackets) using the survey strategy (50, 100'). 





pi 


P2 


P3 


Ae, 




0.702 


-0.474 


-0.530 


0.144 




(0.731)+ 


(-0.435)+ 


(-0.524) + 


(0.164) + 




(0.665)" 


(-0.515)" 


(-0.540)" 


(0.213)" 


r 


0.473 


0.868 


-0.150 


0.043 




(0.412)+ 


(0.895)+ 


(-0.167)+ 


(0.056)+ 




(0.536)" 


(0.833)" 


(-0.134)" 


(0.057)" 




0.531 


-0.145 


0.834 


0.226 




(0.542)+ 


(-0.093)+ 


(0.834)+ 


(0.261)+ 




(0.519)" 


(-0.200)" 


(0.830)" 


(0.327)" 


^-1/2 


0.004 


0.017 


0.271 






(0.005)+ 


(0.020)+ 


(0.313)+ 






(0.006)" 


(0.027)" 


(0.394)" 





certainty (see Tables ^ and 01 , and therefore dominates 
the error on all cosmological parameters. For and erg, 
the MVBs can be approximated using this principal com- 
ponent alone, by crslpjaj sa A8j where j denotes the j^^ 
parameter. This corresponds to considering only the last 
(i = 3) term on the right-hand side of eq. (|SJ). Although 
r is also strongly influenced by P2, this approximation is 
still adequate for this parameter since (T2 is smaller than 
(73 by an order of magnitude. For example, from Tabled 
we infer aslpj^l = 0.163, 0.038, 0.249 for j = 1, 2,3, corre- 
sponding to fim, r, CTg, respectively. These approximations 
are within 10% of the MVB A8j for all three parameters. 

The eigenvectors of the Fisher matrix from the com- 
bined 2PCF ^tot are dominated by - The and (Mfp)- 
eigenvectors show similarities, which can be seen by com- 
paring P3 and Pi for the (f^m, F, erg, JIa) cases (Tables 13 
and IHJi. Accordingly, the largest contribution to ^3(^4) 
comes from o'g(riA), respectively, for both (Affp) and 
In the case of ^tot and the opposite is true. 

For other combinations of parameters, the largest vari- 
ance still makes up more than 80% of the total uncertain- 
ties and also dominates the MVBs. Therefore, an optimi- 
sation scheme should try to maximise the largest variance 
or smallest eigenvalue of the Fisher matrix. This will be 
presented in the next section. 




60 80 100 120 140 60 80 100 120 140 




60 80 100 120 140 



Fig. 3. The variance ai ~ X~ ' of the three eigenvectors 
of the Fisher matrix F, corresponding to (r^m, T, cs) and a 
flat Universe, for (Mfp). Various surveys (TV, R) fSect. l3.1|l 
are compared. The results for the 300 • 13'^ survey is writ- 
ten in brackets in each panel and, while within the range 
of the plot, marked with an arrow. The smooth curves are 
second-order polynomial fits to the data points. 

3.4. Optimisation of the patch radius 

Following KS04, we try to optimise our survey strategies 
by varying the number of lines of sight per patch, N, and 
the patch radius, R, while keeping the total area constant. 
Instead of focusing on the MVB for individual parameters 
as in that previous study, we consider here the eigenval- 
ues Xi of the Fisher matrix. In particular, we concentrate 
on A„, the eigenvalue of the worst constrained principal 
component since this dominates the MVB (see previous 
section). This is sufficient if only three parameters are to 
be estimated from the data, since the variance of the first 
eigenvector by far dominates the others. With four or more 
parameters, however, the second principal component has 
to be included in the optimisation procedure. 

The importance of the principal component corre- 
sponding to the smallest eigenvalue of F on the MVB can 
clearly be seen in Figs. 9-11 of KS04, where the MVB is 
plotted for different N and R. The curves are very simi- 
lar for degenerate parameters, such as (rJm, erg) or (F, rig), 
since these pairs depend on the eigenvectors of F in a 
similar way. Moreover, in the case of a highly dominant 
eigenvector as for the combination (fim, erg, F), all three 
MVBs show the same behaviour, since they all are domi- 
nated by this one principal component. 

In Figs.OHHlwe show the variance ai = 1/A,- corre- 
sponding to the i^^ eigenvector of the Fisher matrix. The 
cosmological parameter combinations are the same than 
in the previous section. We comment on the dependence 
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Fig. 4. The variance ai — \ of the four eigenvectors 
of the Fisher matrix F, corresponding to {flm, T, erg, ^^a), 
for (Mfp). See Fig.|3|for more details. 
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Fig. 6. The variance ai — ' of the three eigenvectors 
of the Fisher matrix F, corresponding to (f^m, T, as) and 
a flat Universe, for ^tot- See Fig.|31for more details. 
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Fig. 5. The variance ai — of the four eigenvectors 

of the Fisher matrix F, corresponding to (f2mj T, erg, Tig) 
and a flat Universe, for (Mfp). See Fig.^for more details. 
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of the Fisher matrix F, corresponding to (rJm, T, o's,ns 
and a flat Universe, for ^tot- See Fig. Ofor more details. 



on the patch radius R and number of lines of sight per 
patch N in the following sections. 

3.4.1. Aperture mass dispersion 

In the majority of the cases, the variance of the dom- 
inating, worst constrained parameter combination shows 
a minimum for some Rq within the probed range of patch 



radii. This confirms the result of KS04, where the MVB 
(which is dominated by tT„) also showed a minimum. The 
optimal radius Rq decreases towards smaller N, thus, 
compact configurations yield better results than sparse 
patches. This is also reflected in the fact that patches with 
high N are preferred over those with small N. For the 
other, less dominant eigenvectors, the corresponding vari- 
ance is (in most cases) a monotonic increasing function 
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Fig. 8. The variance cfi — X~ of the four eigenvectors 
of the Fisher matrix F, corresponding to {Qm, T, erg; ^a), 
for ^tot- See Fig. |21for more details. 



of R. All this implies that in order to obtain constraints 
on parameter combinations, compact, densely sampled 
patches will give the best results. 

3.4.2. Two-point correlation function 

The radius Rq where the variance cr„ of the dominant 
eigenvector takes a minimum is reached for larger R than 
for (M^p). In some cases, Rq exceeds 150 arc minutes, 
which is the maximum of the probed range of radii. In 
most of the time, the variances cr^ are decreasing functions 
of R in contrast to the aperture mass dispersion. These re- 
sults show that a small cosmic variance is more important 
for ^tot than a rigorous sampling of intermediate scales. 
Strategies with sparse patches probing a large number of 
independent regions on the sky and, at the same time, 
capturing shear information on large angular scales will 
provide the best constraints on cosmological parameters. 

We compare the individual contributions of ^+ and ^_ 
to the variances tXi of the joint estimator ^tot- The latter is 
dominated by ^+ and the dependence on the survey strat- 
egy is very similar for both. The contribution from 
however, resembles the one corresponding to the aperture 
mass dispersion, since both statistics sample the conver- 
gence power spectrum in a similar way. In particular, ^_ 
and (Mfp) do not probe large scales in contrast to ^+ (see 
eq. The corresponding filter functions for the former 
two statistics, J4(?y) and [-J i{r]) / rj'^]^ , respectively, are both 
proportional to rj* for small ry, suppressing large scales. 
On the contrary, for the latter statistics, Jq is constant 
for small arguments. A survey which covers large angular 
scales at the expense of a dense sampling of small and 



medium scales will therefore not be optimal for ^_ and 
(A/fp), but efficient for 

3.4.3. Comparison with the uncorrelated 
lines-of-sight-survey 

The 300 • 13' ■^-survey of uncorrelated lines of sight shows 
larger values for ai, i = 2 . . .n than the patch surveys, 
corresponding to poorer constraints on cosmological pa- 
rameters. However, the smallest variance cri is reached 
asymptotically for patch strategies with large R and N — 
10, 20, 30. These surveys consisting of a large number of 
sparse patches corresponding to a small cosmic variance 
are those which are most similar to the uncorrelated lines- 
of-sight-survey. Since for the 2PCF a small cosmic variance 
is crucial, the patch strategies show not much improve- 
ment in (Ti over the 300 ■ 13' ^-survey (unless both ilm and 
r^A are free parameters, see Fig. (S)). This is in contrast to 
the aperture mass statistics, where the improvement in ai 
is more than a factor of two. 

With decreasing j, the variance Ui of the eigenvector pi 
becomes less sensitive on the survey geometry. The mea- 
surement of the best constrained combinations of param- 
eters can therefore not efficiently be improved and even 
the 300 • 13' ^-survey will yield good results. 

The 300 • 13' ^-survey can compete with a patch strat- 
egy regarding the best constrained eigenvectors which con- 
tribute least to the parameter uncertainties. For the domi- 
nant parameter combinations, the patch strategies are su- 
perior and yield much better constraints on cosmological 
parameters. 

3.5. Inclusion of additional parameters 

To include more cosmological parameters in our analy- 
sis, we consider a survey of 141 square degree area. This 
corresponds to an enlargement of the (30, 100') strategy, 
consisting of 100 instead of 10 independent patches. The 
eigenvectors and corresponding eigenvalues are shown in 
Table □ 

The first, best determined eigenvector is orthogonal 
to the f2in-CT8 degeneracy direction, as in Sect. 13.31 The 
second best direction is orthogonal to the prominent T-Ug 
degeneracy direction, but with a strong contribution from 
rim ■ The worst constrained eigenvector is almost solely 
dependent on the redshift parameters zq and /3. 

3.6. Local degeneracy directions 

From the Fisher matrix, we quantify the local direction 
of degeneracy between the parameter pairs (Q^mO's) and 
(r, rig). For each pair, we marginalise over the remaining 
parameters out of {Q^, as,T^ Ug) for a flat Universe, and 
find 

as ^0.56 and UsT^-^ = 0.6, (9) 

assuming a power-law dependence between parameter 
pairs which is usually found in likelihood analysis. These 
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Pi 


P2 


Pi 


P4 


P5 


P6 


P7 






0.639 


-0.563 


-0.016 


0.260 


0.002 


-0.454 


0.014 


0.622 


r 


0.392 


0.691 


-0.404 


0.219 


-0.360 


-0.168 


0.002 


0.233 




0.503 


-0.049 


-0.359 


-0.556 


0.301 


0.463 


-0.034 


0.669 


20 


0.311 


0.088 


0.580 


-0.074 


-0.317 


0.283 


0.611 


4.099 


ris 


0.176 


0.335 


0.266 


0.437 


0.767 


0.078 


0.025 


0.218 


P 


-0.236 


-0.071 


-0.469 


0.011 


0.236 


-0.196 


0.791 


5.292 




-0.038 


-0.279 


-0.285 


0.615 


-0.188 


0.653 


0.004 


0.885 




0.001 


0.004 


0.033 


0.063 


0.117 


1.354 


6.685 





numerical coefficients are basicahy the same for ^tot 
and (M^p). The resulting degeneracy directions are in 
agreement with Kilbinger & Schneider (2005), even 
though a non-Gaussian shear field was used to calculate 
the covariance in this previous work. Marginalisation 
over the hidden parameters increases the volume of the 
error ellipsoid and also alters the orientation of its axes 
in parameter space. However, the general degeneracy 
direction is similar for various marginalised parameters, 
see Figs. El and ^1 We conclude that the directions of 
near-degeneracy between the considered parameter pairs 
are robust against the inclusion of non-Gaussianity of 
the shear field, but less stable against the addition of 
cosmological parameters. In Fig.lHlwe show the Icr-ellipses 
for the parameter pair (rim, F) using all four estimators 
^-1-, ftot and (Afj^p). Slight misalignments in the 
degeneracy directions from ^+ compared with ^_ leads to 
improved constraints from the combined estimator ^tot- 



3. 7. Effect of the survey depth on the principal 
components 

For surveys with fixed A'^ = 60 and varying R, we calcu- 
late the Fisher matrix for three different source redshift 
distribution functions, parametrised by zq = 0.8, 1.0 and 
1.2, respectively. The variance ai of the principal com- 
ponents of the Fisher matrix corresponding to the three 
parameters (rimjCs, F) is shown in Fig. ^2 The general 
characteristics of Aj as a function of R does not change 
with the survey depth and simply gets scaled. The vari- 
ance corresponding to the first two eigenvectors decrease 
with increasing survey depth, corresponding to a smaller 
error on the parameter measurement. However and un- 
expectedly, the principal component corresponding to the 
largest variance 173 is best determined for the most shallow 
survey with zq = 0.8. This principal component takes little 
contribution from F and points in the Q^n-'^s degeneracy 
direction. We repeat the PCA without the shape param- 
eter and found a similar result. Also when marginalising 
over additional parameters, the largest eigenvalue of the 
Fisher matrix is smallest for zq = 0.8. 

This unexpected effect is a local one, present in the 
Fisher matrix only. In order to see the global behaviour, we 




Fig. 9. la error ellipses in the f7m-F plane. From 
small to large ellipses, successive marginalisation over 
(78 5 ^s, f^A, Pi Zq (except the ones that are plotted) was per- 
formed. The four panels correspond to the four estimators 
^tot and (Afap) as indicated. 



calculate the likelihood in the rim-cs plane and find that 
the confidence levels get tighter with increasing source 
redshift. The near-degeneracy is less pronounced if the 
survey is deeper and more information about the large- 
scale structure is collected in the shear signal. Moreover, 
the curvature of confidence levels increases with increas- 
ing Zq which leads to a larger bending of the curves of 
constant likelihood. The shear correlation from different 
source galaxy redshift distributions allows one to constrain 
slightly different regions in parameter space. This fact is 
made use of in shear tomography to lift the parameter 
near-degeneracies (Hu 1999). 
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Fig. 10. 1(7 error ellipses for and erg, rig, fl^ and T, re- 
spectively, in the four panels. From small to large ellipses, 
successive marginalisation over the remaining parameters 
out of (r, (jg, rig, i^A, ^o) was performed. The estimator 
is (M^p) and the survey strategy the 141 square degree 
survey. 



3.8. Effect of the source galaxy density on the 
principal components 

For the same survey types {N = 60) and the parame- 
ters J7m,f8 and F as in the previous section, we com- 
pute the Fisher matrix for three different surface densities 
of background galaxies, rigai — 20, 30 and 40 per square 
arc minute. The variance ct^ of the three principal com- 
ponents of the Fisher matrix is plotted as a function of 
patch radius R in Fig.^l corresponding to measurements 
of the 2PCF and the aperture mass dispersion, respec- 
tively. As expected, ai decreases with increasing rigai since 
more background galaxies provide a better sampling of the 
shear field. The minimum variance for era at Ro ~ 130 in 
the case of (Affp) does not change with noise level. 

In this work, we treated the survey depth, 
parametrised by zq , and the number density of background 
galaxies rtgai as independent survey characteristics. For a 
realistic survey however, this assumption is not justified 
in general since an increase in depth implies a higher 
galaxy number density. 



3.9. Principal components and the error ellipsoid 
volume 

We investigate the geometric mean p — (Ai . . . A„)^/" of 
the n eigenvalues of the Fisher matrix F. The equiva- 
lent radius of the Fisher error ellipsoid, which is the ra- 
dius of an n-sphere with the same volume, is p~^/'^ = 
(detF)-i/2n^ 
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Fig. 11. The variance at = of the three eigenvectors 

of the Fisher matrix corresponding to (fJnijfSjr) and a 
flat Universe, using the 2PCF (top) and the aperture mass 
dispersion (lower panels). The curves show different survey 
depth where solid, dashed and dotted lines correspond to 
Zq = 0.8, 1.0 and 1.2, respectively. 

In Figs. El and E| we plot the effective error ellipse 
radius p~^/^ as a function of R for survey geometries 
{N, R) using different parameter combinations and esti- 
mators. Since p~^l'^ is dominated by the largest variance 
(T„ (smallest eigenvalue of F), it shows a similar behaviour 
(see Sect. l3.4|l . However, the minimum seems to be reached 
at smaller i?o an therefore the optimal survey radius is 
smaller than the previous sections implied. The 300 -13'^- 
survey (see Fig. E| and Table IHl yields comparable re- 
sults than for the patch geometries if ^+ is used. ^_ and 
(M^p) (not shown), however, strongly suffer from the lack 
of scales larger than 20 arc minutes which results in a 
much larger error ellipsoid. 



Table 8. The effective radius of the Fisher error ellipse 
p~^l'^ for the 300 • 13'^-survey, corresponding to two com- 
binations of cosmological parameters. The values of p~^l'^ 
for various patch strategies (A^, i?) are displayed in Fig. 1141 
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Fig. 12. The variance ai = Aj of the three eigenvectors 
for the Fisher matrix corresponding to (r2m,cr8,r) and a 
flat Universe, using the 2PCF (top) and the aperture mass 
dispersion (fower panels) . The curves show different survey 
depth where solid, dashed and dotted lines correspond to 
fT-gai — 20, 30 and 40 per square arc minute, respectively. 



In Fig. El the results for the 2PCFs £,+ and are 
compared. As already mentioned before, ^tot is dominated 
by ^+ and shows a similar behaviour. The shape of the ^- 
curves on the other hand resembles the one for (M^ 



3.10. Principal components of ttie scaled Fisher matrix 

The scaled Fisher matrix F is defined as 



F 



(10) 



31 



Its inverse F"^ represents the covariance of the nor- 
malised (unit variance), correlated Gaussian variables 
constructed from the original cosmological parameters. 



9j = ^il ji- The scaled parameters can be compared 
with each other in a straightforward way, and their cor- 
relation gets more underlined. The variation with patch 
radius R of the eigenvalues of the scaled Fisher matrix 
matches well with their unsealed counterparts. In Fig. 1151 

2/2 

we plot the variance of the dominant eigenvector A„ (i?) 
for geometries with = 10 and A^ = 60, respectively, and 
the cosmological parameters (f2m, T, cs, ^k). 
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Fig. 13. The effective radius of the Fisher error el- 

lipsoid as a function of survey radius i?. The left panels 
correspond to (Af^p), the right panel to ^tot. In the up- 
per two panels, the parameter combination is (rim, T, erg), 
in the lower panels it is (r2,n, F, (Tg, JIa)- The numbers in 
parenthesis correspond to p^^l"^ from the 300- 13' ^-survey. 
The definition of p is given in the panel as a reminder. 
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Fig. 14. The effective radius p^^l"^ of the Fisher error 
ellipsoid corresponding to (rii„,r, cg) (left panel) and 
(rim, r, cg, riA) (right panel). The three estimators 
^+ and ^tot are compared using for surveys with sparse 
(iV = 10, solid lines) and dense (A^ — 60, dash-dotted 
lines) geometries. 



3.11. Principal components and the use of priors 

Additional priors on parameters modify the original Fisher 
matrix F to F — F + C~^. For example, a Gaussian 
prior for the i'^ parameter with variance Si corresponds 
to C^^ = Sij s,~^. Priors lower the eigenvalues. Since the 
reduction is the larger the smaller the eigenvalue is, the 
dominant eigenvector with the smallest eigenvalue is af- 
fected most by priors. The effect of a prior on the vari- 
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Fig. 15. The variance — of the dominant eigen- 

value of the scaled Fisher matrix F for the parameters 
(rim, r, fS) ^a)j as a fmiction of the patch radius R. Solid 
lines correspond to surveys with N — 10, dash-dotted 
curves are for N = 60. 
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Fig. 17. The variance of the dominant eigenvector CT4 = 

— 1/2 

A4 of the scaled Fisher matrix F for the parameters 
(fijn, F, as, ng), as a function of the patch radius R. Thick 
curves are for a flat Universe, thin curves for constant fl^. 
Solid lines correspond to = 10, dash-dotted lines to 
N = 60. 
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Fig. 16. The variance of the dominant eigenvector 0-3 is 
plotted for patch geometries with = 10 (solid lines) 
and = 60 (dash-dotted), for (Mfp) (left panel) and ^tot 
(right panel). The three cosmological parameters under 
consideration are fim, F and as- The upper curves cor- 
respond to the case without prior. For the lower pair of 
curves, and additional prior on F of s(F) = 0.3 (0.05) is 
assumed for (Mfp) (^tot), respectively. The values in the 
parenthesis denote cr„ the 300 ■ 13' ^-survey. 



ance of this eigenvector is shown in the upper right panel 
of Fig.^and in Fig.^] In Fig.^ the priors s{pi) are 0.003 
for riniiCTg and F, 0.03 for Oa, 0.1 for rig and Zq and 1 for 
Pp. For the three curves, priors are added for parameters 
as indicated in the panel. The more priors are included, 
the more is the large, dominant variance affected, and the 
curve ai as a function of i becomes less steep. 

The addition of a prior also flattens the curve cr„(i?) 
as a function of patch radius R and the difference between 
the surveys becomes less pronounced, making the optimi- 
sation less effective, see Fig. ^| From Fig. El one infers 
that the prior of a flat Universe (f^A = 1 — 51i„) corre- 
sponds to a higher variance (T„ than a constant Q\ — 0.7 
prior. 



4. Summary and discussion 

We study the design of cosmic shear survey with the means 
of a principal component analysis of the Fisher informa- 
tion matrix. Although only a local approximation of the 
likelihood, the Fisher matrix provides a simple and cost- 
effective method to predict cosmological parameter con- 
straints and to investigate degeneracies between parame- 
ters. We compare realistic survey designs with non-trivial 
topology. In most of the previous analyses (an exception 
being KS04) either simple monolithic survey topologies 
were assumed or the correlation of power for different 
Fourier modes was neglected. In this work, we use the 
full covariance matrix of second-order shear statistics cor- 
responding to a complicated, non-trivial distribution of 
lines of sight. We assume the shear field to be Gaussian 
which leads to an underestimation of the covariance on 
intermediate angular scales between ~ 1 and 10 arc min- 
utes. 

The eigenvectors or principal components of the Fisher 
matrix determine linear, orthogonal combinations of cos- 
mological parameters. We find that these vectors are only 
little affected by the survey characteristics. The same is 
true for the directions of near-degeneracy between pa- 
rameters. However, the eigenvalues change significantly 
with survey geometry. The eigenvector of the Fisher ma- 
trix which dominates the errors of the cosmological pa- 
rameters has the smallest eigenvalue corresponding to the 
largest variance. To maximise this eigenvalue, which is to 
minimise the error on this eigenvector, the shear correla- 
tion on large angular scales up to several degree has to be 
measured. The dependence on the sampling of shear cor- 
relation is different for the correlation function ^tot and 
the aperture mass dispersion (Mfp). For ^tot, a survey 
consisting of sparse patches with a small cosmic variance 
yields the smallest errors. Using (Mfp), dense and slightly 
smaller patches are preferred. 

We consider various combinations of the cosmological 
parameters rim, T, as, rig, a and the source galaxy red- 
shift characteristics zo and /3 (see Sect. I2.l(l . If not more 
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than three parameters are to be determined from our 
survey (the other parameters being fixed), the dominant 
eigenvector of the Fisher matrix comprises more than 99 % 
of the parameter errors. For more than three free param- 
eters, the first two principal components are responsible 
for most of the parameter errors. 

Changing the number density of source galaxies, the 
depth of the survey, the intrinsic ellipticity dispersion 
or the total survey area simply causes a scaling of the 
eigenvalues. The optimal survey setting is not affected 
by changes in those parameters. Furthermore, introduc- 
ing external priors on some parameters does not result 
in a different optimal setting. However, this optimum be- 
comes less pronounced so that optimisation of the survey 
design will be less important. If a flat Universe is taken as 
a prior, the dependence on the survey geometry is much 
weaker than for the prior Q\ = const. 

The dominant principal component shows an unex- 
pected behaviour when the survey depth zq is varied. 
The error on the corresponding parameter combination 
increases for increasing zq in the range of zq = 0.8 ... 1.2. 
However, this effect is only local and disappears when con- 
sidering the likelihood function, which shows strong non- 
Gaussian features not present in the Fisher matrix. 

As survey sizes increase it will be desirable to estimate 
an increasing number of parameters independently from 
weak lensing surveys alone such as the equation of state 
of dark energy or its variation with redshift (Hu 1999; 
Heavens 2003). However, parameter near-degeneracies in- 
herent in weak lensing observables make additional data 
from different cosmology experiments very valuable in 
breaking such degeneracies. On the other hand, priors 
from CMB or SN observations can be used in the future 
to optimise the design of smaller surveys by using princi- 
pal components of the joint Fisher matrix (e.g. Crittenden 
& Pogosian 2005 provide results of joint PCA of different 
data sets in a similar context). 

We conclude that existing tools such as generalised 
eigenmode analyses (Karhunen 1947; Loeve 1948; Vogeley 
& Szalay 1996; Matsubara et al. 2000; Kilbinger & Munshi 
2005) and principal component analyses which we have 
studied here will be very useful in optimising future weak 
lensing surveys. 
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